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ABSTRACT 

Three-minute oscillations over sunspot's umbra in AR 11131 were observed 
simultaneously in UV/EUV emission by SDO/AIA and in radio emission by 
Nobeyama Radioheliograph (NoRH). We use 24-hours series of SDO and 8-hours 
series of NoRH observations to study spectral, spatial and temporal variations 
of pulsations in the 5-9 mHz frequency range at different layers of the solar 
atmosphere. High spatial and temporal resolution of SDO/AIA in combination 
with long-duration observations allowed us to trace the variations of the cut-off 
frequency and spectrum of oscillations across the umbra. We found that higher 
frequency oscillations are more pronounced closer to the umbra's center, while 
the lower frequencies concentrate to the peripheral parts. We interpreted this 
discovery as a manifestation of variation of the magnetic field inclination across 
the umbra at the level of temperature-minimum. Possible implications of this 
interpretation for the diagnostics of sunspot atmospheres is discussed. 

Subject headings: sunspots - Sun: oscillations - Sun: UV radiation - Sun: radio 
radiation - Waves 
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Introduction 



Despite the fact that three-minute sunspot osciUations have been observed for 
many decades and enormous theoretical process has been made in understanding this 



phen omenon (for reviews see 



Bogdan fc Judge 



2006 



Thomas fc Weiss 



2008 



Khomenko 



20091 ) . there are stiU many questions concerning sunspot wave physics. It is generaUy 
accepted, that they are slow magnet oacoustic waves moving with the local sound speed 
along magnetic field lines from the photosphere through the chromosphere into the corona. 
However, it is still unclear, how high in the corona they can propagate? What determines 
their properties in the corona? What is the spatial variation of the cut-ofl[ frequency 
across the umbra? Are they directly linked with the propagating longitu dinal waves 



ubiquitously observed in th e EUV fans of coronal active regions (e.g. see 



De Moortel 



200^; 



Verwichte et al. 



2010l . for recent reviews). What is the nature of the recently 



revealed relationship of t 



releases (Sych et al. 



le three-min oscillations and similar periodicities in fiaring energy 



20091 )? Understanding these issues provides an opportunity for plasma 



diagnostic at different layers of solar atmosphere since these waves conve y the infor n iation 



about the properties of their formation and propagation mediums. Thus, 



Shibasakil fcoOlh 



proposed that umbral oscillations seen in microwaves can be used as a diagnostic of the 



temperature-minimum region. Further 



Solankil ( 120031 ) pointed out that observations with 



higher spatial resolution are needed to isolate the umbra and to improve the diagnostics. 

Spatial distribution of 3-min oscillation frequency and power across a sunspot 
at different heights has been intensively studied with the use of data obtained with 
ground-based and space-borne instruments (e.g. onboard SOHO, TRACE and Hinode). 
In particular, SOHO/MDI observations showed that at the photosphere the oscillation 
power is e nhanced in the Dop pler velocity and line-depth signals, but not in the continuum 



intensity ( jjain fc Haber 



2OO2I ) . Enhancement of the spectral power was found in the 
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regions of reasonably strong magnetic field, but showed evident suppression in regions 
of strongest fields, in the center of the sunspot. Observations in demonstrated that 
standing oscillations with the frequencies around 6.5 mHz were dominant in the umbral 
and inner penumbral regions, together with the presence of running penumbral waves in 
the band of 3 mHz. However th ese waves were not detected at the photospheric heights 



( jChristopoulou et al. 



2000 



200 ll ). It was suggested that both 3-min umbral oscillations 
and 5-min running penumbral waves were in fact different manifestations of the same 
wave, generated in the photosphere, and propaRatinR alouR di f ferent magnetic field li nes 



with the different inc 



Christopoulou et al. 



inati on angle ( Ivan der Voort et al. 



2003 



Bloomfield et al. 



2003). 



( 120031 ) demonstrated that oscillations with different frequencies in 



the 3-min band were localized within different parts of the umbral chromosphere. Then 



Tziotziou et al. 



( 120061 ) found that in the chromosphere, the frequency of 3-min oscillations 
decreases from the umbra outwards, with the abrupt change at the umbra-penumbra 
boundary. A small "node" at the center of the umbra, wh ere the oscillation s with the 



frequencies above 4 mHz were suppressed, was reported by 



Nagashima et al. 



(120071). A 



similar eff ect of "calmest umbral position" was found in the Doppler velocity maps in Ca H 



8542 A by 



Tziotziou et al. 



( 120071 ). Nature of this "black spot" in the spatial distribution of 



the chromospheric 3-min oscillations remains unrevealed. 

In the transition region the 3-min oscillations are reported to fill in all the sunspot 



umbra. In the corona t hey concentrate in sma l. 



of bright coronal loops ( IBrynildsen et al. 



2000 



er regions t h at coincide with t 



20021 . 



2004h . 



leendpoints 



De Moortel et al. 



( I2OO2I ) used 



TRACE observations of EUV intensity waves in 38 coronal loops and found that in the 
loops situated above sunspots the waves have periods of about 3 minutes. This suggests 
that the underlying oscillations can propagate through the transition region into the corona. 



Hence, it is still not clear how the oscillation power and frequency in the 3-min band 
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are distributed within the umbra of the sunspot at different altitudes. This knowledge is 
important for revealing the nature of 3-min oscillations, and for diagnostic of the solar 
atmosphere above sunspots. 



The data from the Atmospheric 



Solar Dynamics Observatory (SDO; 



Ima^inR Assembly (AIA; iTitle et al.ll2006l ) on the 



Schwer et al. 



20021 ). with their high resolution (0.6 



arcsec per pixels) and high cadence (up to 12 s), offer a new opportunity for such a 
study. The instrument observes solar plasma from photospheric to coronal temperatures, 
taking full-disk images in a variety of extreme ultraviolet (EUV), ultraviolet (UV) and 
visible light wavelength bands. The other advantage of SDO data is 24-hour continuous 
observations. Three-min oscillations have very high quality, they can be observed for many 
days. Therefore, long-duration observations allow us to study this phenomenon with very 
high spectral resolution, which has not been done before. 

In this paper we analyze simultaneous observations of the sunspot-associated 
oscillations by two instruments, AIA and the Nobeyama Radioheliograph, in order to obtain 
the information on spectral, spatial and temporal structures of oscillations above the same 
sunspot. We restrict our analysis to the frequency range 5-9 mHz and horizontal extent to 
the umbral area. 



2. Observations 

Active region (AR) 11131 was available for observations in the northern hemisphere 
near to 30 degrees latitude and crossed the central meridian 2010 December 8. We 
selected this date for our analysis. Figure [I](a) presents line-of-sight magnetogram of the 
sunspot region from Helioseismic and Magnetic Imager (HMI) on SDO. AR consisted 
of a well-developed sunspot of southern polarity with the magnetic field strong enough. 
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B = —2300 G to produce significant and highly polarized emission at 17 GHz (A = 1.76 
cm). The sunspot was surrounded by a number of small magnetic fragments of northern 
polarity. 

The emission mechanis m of radio sources a ssociated with sunspots is known to be 



gyroresonance emission (e.g. 



Zheleznyakov 



19621 ). At 17 GHz, the third harmo nic layer i s 



the m ost probable emission region, corresponding to a 2000 G isogauss layer ( jShibasaki 



19M) 



The radio map of the sunspot that we studied obtained by Nobeyama Radioheliograph 
(NoRH) is shown in Fig. [I](b). Color indicates brightness temperature at 1.76 cm (Stokes 
I). NoRH measures both right circular polarization and left circular polarization. Using 
these measurements we synthesized a set of NoRH images in I and V Stokes parameters 
with time interval of 10 s. The NoRH beam width (FWHM) is indicated by the orange 
ellipse at top-right. Due to low elevation of the Sun during the winter time, the beam size 
was roughly 17^' x 30^', comparable with the radio source, and therefore, determined an 
elliptical shape of the latter. 

In Figure [I](c-h) we show intensity images of the sunspot in a sequence of increasing 
temperatures from six AIA bandpasses: 1700 A, 1600 A, 304 A, 171 A, 193 A, and 211 A. 
The spatial resolution of data is 0.6 arcsec per pixel. 

Atmospheric Imaging Assembly observes UV continuum in 1700 A bandpass with the 
characteristic formation temperature about 5000 K. This bandpass corresponds to the lowest 



formation layer located close, but sti. 
to the sunspot atmospheric model by 



1 above the t e mper ature-minimum region. According 



Maltbv et al 



( 119861 ) , the temperature-minimum region 



is located at about 500 km above t^qq = 1, corresponding to the upper photosphere/lower 
chromosphere. The AIA 1600 A bandpass has a broad response to temperatures ranging 
from the temperature-minimum to chromospheric line C IV emissions (4 — 10 x 10^ K) 
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in active regions, but the dominant contributor is the UV continuum above the region of 
temperature-minimum. The 304 A channel is dominated by the two He II hues with average 
formation temperature T = 5 x 10^ K, corresponding to the transition region. EUV hght 
from the Fe IX emission hne in the 171 A bandpass primarily coming from upper transition 
region/quiet corona material with average formation temperature T = 6.5xlO^K. The 
193 A channel is expected to be dominated by Fe XII line (average T = 1.5 x 10^ K) and 
211 A channel is exp ected to observe Fe XIV line (T = 2 x 10^ K) for AR observations 



to'Dwyer et al 



2010|). 



We did not use the other SDO/AIA bandpasses for the following reasons: visible light 
continuum observations at 4500 A are available only with 3600 s time cadence; observations 
in coronal channels 94 A and 131 A had not sufficient quality for this sunspot. 



3. Data analysis 



3.1. Temporal variations of UV/EUV and microwave emission 



It is known that changes in the UV/EUV emission with time must be due to changes 
in the temperature and density in the emission forming region. Plasma temperature and 
density variations, in turn, are thought to be evidence some dynamical processes, e.g. a 
compressible wave passage through the corresponding region. 



Different interpretations for the radio brightness oscillation were proposed. 



as beiuR caused by the variations in magnetic field strength ( IGelfreikh et al. 



'or ex ample, 



19991 ). 



Shibasakil (120011 ). however, showed that they can be explained consistently with the EUV 



observation in terms of the density and temperature fiuctuations due to upward-traveling 
acoustic waves through the third harmonic gyroresonance layer (2000 G for observations at 
1.76 cm). 
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To study temporal variations of UV/EUV emission above the center of the umbra, 
where the magnetic field is presumably vertical, we analyzed de-rotated 24-hour series 
of intensity obtained from one pixel, corresponding to the umbra's center and shown by 
white cross on the 1700 A image of Fig. [11(c). The neighbouring pixels show a very similar 
time variability, hence the data are not affected by the possible variation in pointing. For 
microwave emission we picked up the value of the maximum brightness temperature, TBmax) 
located in the center of the elliptical source and, therefore, corresponding to the sunspot 
umbra center for every synthesized map. 

Figure [2](a) shows time profiles, /(t), of the UV/EUV intensity (in arbitrary units) 
and microwave peak brightness temperature (in K) of 22 minutes duration for clarity. 1600 
A intensity is multiplied by a factor of 3 and NoRH 1.76 cm brightness temperature by a 
factor (4 x 10"^). Time cadence is 24 s for 1700 A and 1600 A, 10 s for 1.76 cm and 12 s 
for other channels. Well-pronounced variations of the emission intensity, with the periods 
of about two-three minutes, are clearly seen in all channels. 

The time profiles of the modulation depth are shown in Fig. El^b) for 1700 A, 1600 A, 
304 A, and 1.76 cm and in Fig. El^c) for 171 A, 193 A, and 211 A usine: the same colors as 
on the top panel. This quantity is calculated as 

A/// = {I{t) - /o)//o, (1) 

where I{t) is time series of intensity of UV/EUV or microwave brightness temperature. The 
slowly varying mean signal /q is obtained by smoothing the signal I{t) over an interval of 14 
minutes. The highest modulation depth is found for 304 A emission, 20% on the average, 
and it decreases down to the photosphere, 10% for 1600 A, 8% for 1700 A, and up to the 
corona with hight, 18% for 171 A, 12% for 193 A, and 8% for 21lA. 

Modulation depth 5% is found for microwave brightness variation. It corresponds to 
the change in brightness temperature by ATg 1.2 x 10"^ K. These values agree with 
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previous results showing 1 - 5% of Tb variation. 



The general trend of the height change in modulation depth of UV/EUV emission also 



agrees with results of previous studies, which have shown that t. 



peaks at transition-region temperatures (e.g., 



O'Shea et all 



le oscillation amplitude 



2002 



Marsh fc Walsh 



20061). 



However, it is rather surprising to find intensity variations with such a high amplitude in 
coronal channels. Previous studies showed umbral oscillations, observed with TRACE 171 A 
and 195 A, with much smaller amplitudes (3-8%) usually guided along the EUV fans of the 
limited transverse size. One possibility is the presence in the coronal EU V channels of AIA 



Brynildsen et al. 



of the emission generated at much lower temperature and hence heights. 
( 120041 ) estimated the contributions from the transition region lines to the TRACE 171 A 
and 195 A channel intensity. They found an average value of this contributions equal to 
17% and finally concluded that the significant part of th e oscillations in the T RACE 171 



Q^Dwyer et al. 



( 120101 ) performed 



A and 195 A channels must indeed occur in the corona, 
similar examination for AIA and found that the dominant contribution to three EUV 
coronal channels under study from the AR plasma comes from the expected ion line formed 
at the hot, coronal temperatures. 



3.2. Cross-correlation analysis 

The other indirect indicator of the brightness fiuctuation origin in each AIA channel 
can be time shifts between corresponding signals. If the pulsations originate from the same 
layer of the solar atmosphere we may expect an in-phase behavior for intensity variations. 

Cross-correlation coefficients as a function of the lag for each analyzed wavelength 
with respect to 304 A intensity variation are shown in Figure [3](a) and corresponding 
delays in Figure [3](b) . We used temporal interpolation of the 304 A time series to match 



-10- 



other channel observing times. The precision hmit of the obtained delays is estimated by 



a = k X where = y^At^/T2 is the standard deviation of the sampling error for the 
uniform time resolution At and k = 2 for 95% confidence from the Student-t distribution. 
These uncertainties ±cr are shown by error bars in Figure [3](b) . 

Time delays increase with the expected formation height of different channels, except 
211 A channel. Equal delay for 171 A and 211 A (12^ ± 7^) supposes that possibly 211 
A intensity has a contribution of cooler temperature lines than expected, or the cooler 
and hotter plasma structures simply co-exist in the corona at the same height. At the 
same time, increasing delays found for other wavelength suggest that observed intensity 
variations in different channels originate from progressively higher altitudes and result from 
waves propagating upwards along magnetic field lines. The cross-correlation coefficient for 
1700 A bandpass is poor, probably due to very weak emission variation inside the umbra 
region and nonuniform temporal resolution of downloaded data, and is not shown here. 

If we estimate the expected height separation AH between the two consecutive AIA 
wavelengths based on the mean formation temperatures mentioned in Section O we can 
use obtained delays to make rough estimates of the average velocity of wave propagation 
between the two heights. Assuming 1600 A formation height around temperature-minimum 
(about 500 km) and 304 A in the transition region (about 2200 km), AHi ^ 1700 km 
gives vertical velocity Vi ^ 70 km s~^. This value exceeds the speed of sound in the 



chrom osphere, which lies between 10 and 30 km s ^. However, 



Abramov-Maximov et al. 



(I2OIII ) found a close value of speed (60 km s ^) based on cross-correlation between optical 



and radio oscillation. 

The height separation AH2 ^ 1000 km between the heights of the 304 A and 171 A 
emission layers gives velocity V2 ^ 83 km s~^. This value corresponds to the calculated 
sound speed with T = 2 x 10^ K, that is an average temperature between the two formation 
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layers. 

It is rather difficult to define the height of 193 A emission formation over the sunspot 
umbra. Bearing in mind the speed of sound (130 km s~^) calculated for the 171 A emission 
formation temperature (6.5 x 10^ K), the delay 12^ ± 7^ between the 171 A and 193 A 
emissions gives the minimum value of AHs > 1500 km. Therefore, the 193 A emission 
above the umbra corresponds to the corona with the T > 1 MK. 

Cross-correlation analysis gives 10s ± 7^ of the time delay between the 304 A and 
microwave signals, which means that the formation height of radio emission is between 
AIA 304 A and 171 A bandpass. Using V2 ^ 83 km s~^ obtained for this height interval, 
we estimated the altitude of 1.76 cm emission generation, which is approximately 3000 
km. The observed value of peak brightness temperature Tsmax = 2.5 x 10^ K exactly 
corresponds to the estimated altitude. 

We have to note that the time delays presented in Figure [3] are averaged over the 
frequencies and time. Further analysis reveals that these delays are frequency- and 
time-dependent. In Figured] we show the emission modulation depth of signals at 304 A 
(red) and 193 A (black) at the center of the umbra, shifted in time by 5t to match each 
other. Signals are filtered in narrow (1 mHz) frequency bands with the central frequencies 
at 5 mHz, 6 mHz, 7 mHz and 8 mHz (± 0.5 mHz). In this example the 193 A signal is 
delayed against 304 A signal by 40 s at 6 mHz and 7 mHz, and by 24 s at 8 mHz. But at 
5 mHz we found an inverse delay: the 304 A signal is delayed against the 193 A signal. 
Moreover, these phase delays are not only frequency-dependent, but also time-dependent: 
though varying the time shift 6t we could adjust the phases of the signals in one wave train, 
the same phase shift did not do it in the other wave train. The phase shifts were found 
to vary strongly during the 24-hour observation. Often it was not easy to resolve the 27r 
uncertainty. 
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The fact of the frequency dependence of the delays between the signals measured at 
different heights agrees with previous empirical results obtained from the c ross-correlation 



of photosph eric and chromospheric velocity maps filtered in narrow bands (ICenteno et al. 



20061 . 



20091 ) . as well as with the theoretical prediction that oscillation modes in dispersive 
medium propagate with the phase velocity depending on the frequency. The time 
dependence i s also consistent with the previous results, e.g. the frequency drifts in the 



3-min band (jChristopoulou et al. 



20031 : 



Sych et al 



201ll ). Therefore, detailed analysis of 



the phase relations between the perturbations at diflFerent heightsis of obvious interest. 
However, it is beyond the scope of this paper and will be addressed in a follow up study. 



3.3. Spectral analysis 

To study spectral characteristics of intensity variation we used the fast Fourier 
transform (FFT) method. Because of long intervals of uniform data with high cadence, 
the spectral analysis was done with high-frequency resolution, 36 /xHz for NoRH and 12 
//Hz for AIA. Such high spectral resolution in combination with the high spatial and 
temporal resolution of SDO/AIA allowed us to trace the variations of the spectrum of 3-min 
oscillations across the umbra. In particular, we were able to trace the so-called low cut-oflF 
frequency, which is defined here as a lowest frequency in the 4-9 mHz band for which the 
Fourier power is three times higher then the noise level. 

Figure [5] represents the power spectra covering 24 hr of UV/EUV data and 7 hr 
45 min of microwave data. The power is normalized on its maximum in the band 4-10 
mHz. The spectra are shown for the single pixel above the center of the umbra. Most 
of spectra have dominant frequencies between 6 mHz (period 167 s) and 7 mHz (period 
143 s). High-frequency tail extends beyond 8 mHz. All UV/EUV spectra exhibit a sharp 
increase at 5.7 mHz (equivalent period 175 s), which can be considered as the lower cut-off 
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frequency. 



It is interesting that a similar value (5.8 mHz) was recently found for the maximum 



cut-off frequency at the axis of a model sunspot in the 3-D nonlinear MHD simu 



the driver located under the photosphere and emitting 5-min waves (IFelipe et al. 



ation with 



2010h . 



The spectrum of the radio oscillation at 1.76 cm does not show a clear low cut-off. 
However, the microwave signal has a duration of less than 3 times that the UV and EUV 
data and comes from much larger area: the difference in the pixel size of AIA and NoRH 
maps is about 10 times. Nevertheless, similar to all AIA spectra, it has a maximum between 
6 and 7 mHz and the increase in the power near 8.8 mHz. 

In previous studies the value of the cut-off frequency was often defined as the highest 
frequency where the phase difference between two signa ls measured at different heights was 



zero (e.g.. .Centeno et al, 



20061 



2009: 



Felipe et ah .2011 ). As we see from Fig. (U the phase 



difference between the 193 A and 304 A is positive in the 6-8 mHz band, while is negative 
at 5 mHz. This confirms that cut-off frequency is located between 5 and 6 mHz. However, 
the above-mentioned method allows one to determine the cut-off frequency averaged over 
the umbra, while we are interested in the spatial distribution of this value across the umbra. 

Further analysis has shown that spectral characteristics of oscillations vary across the 
umbra. Figure [6] presents the spectra for pixel located further from the center, near the 
umbra/penumbra border, it is shown by white asterisk on Fig. [I](c). It is clearly seen that 
the low frequency cut-off dropped down to 5 mHz and the peak frequency moved down to 
5.5 mHz as compared with the umbra's center. The same tendency is found for all other 
AIA channels. 

To trace the changes in the spectrum of oscillations for the entire sunspot region we 
calculated Fourier power spectra for every spatial pixel of the intensity maps. Figure [71 
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gives the spatial distribution of Fourier power in four narrow (A/ = 0.4 mHz) frequency 
bands with the central frequencies at 5 mHz, 5.5 mHz, 6 mHz, and 7 mHz (±0.2 mHz in 
each band). Individual power spectra from each spatial pixel have been normalized on 
its spectral maximum. Normalization was performed to aid in the comparison of spectral 
profiles between spatial regions that exhibit vastly different intensity amplitudes. Only 
pixels with power larger than three times the noise level are plotted in the power spectrum 
maps. The umbra-penumbra boundary as observed in 1700 A continuum is shown in every 
plot by the contour line, white in panels (a-b) and black in panels (c-f). The external 
penumbra border is also shown in the upper left plot of panel (a) for reference. 

It is clearly seen from Figure [71(a) that Fourier power in all four frequency bands is 
confined inside the umbra and avoid the penumbral region. In the 5 mHz band it also avoids 
the central part of the umbra and forms a ring near to the umbra/penumbra border. In the 
5.5 band this ring extends inside the umbra, becoming wider. We already saw in Figure 
[5] that oscillations at frequencies / < 5.7 mHz are absent in the center of the sunspot. 
Oscillations in 6 and 7 mHz bands cover the whole umbra, but are more pronounced in its 
central part. This tendency is seen through all the AIA bandpasses, see Figure [7l(b-f). At 
the same time the oscillating area slightly broadens with the height. The diameter of the 
oscillating circle is widened by about Ad = V with each subsequent bandpass, expanding 
beyond the umbra boundary at 304 A. In coronal channels 171 A - 211 A pulsations with 
high power outspread to the coronal fan-structures. Higher frequencies, 8 and 9 mHz exhibit 
even stronger concentration toward the center, but lower power (they are not shown here). 

The peculiarities of the spatial distribution of 3-min oscillation frequency found here 
are consistent with previously obtained results. In particular, with the decrease in the 
peak frequency from the umbra outwards wi th the abrupt change i n the umbra-penumbra 



boundary in the chromosphere, reported by 



Tziotziou et al. 



(120061 ). Also, the existence 
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of the region with the su ppressed powe r of 3 -min osciUations at the center of the 



umbra has been found by 



Jain fc Haberl (|2002r) in th e phot ospheric Doppler velocity and 



hne-depth data of SOHO/MDI, by 



m Hz) usinR Hinod e /SO 



Nagashima et al. 



( 120071 ) at lower frequencies (5-5.5 



by 



Tziotziou et al 



observations of low-chromospheric Ca II H line intensity, and 
( I2OO7I ) in their Doppler velocity maps obtained at Ca II 8542 A. Here 
we show that such a region at 5-5.5 mHz exists also at the transition region and coronal 
heights, and that its size is frequency-dependent. 



4. Interpretation of spectral characteristics 

It is interesting that spectral characteristics of 3-min oscillations, namely, the low 
cut-off frequency and its variation across the umbra, are similar through all the temperature 
bandpasses. This finding suggests that the strong magnetic field of a sunspot works as a 
waveguide for the acoustic waves generated already at the photosphere level and eventually 
reaching 1 MK corona. 

To understand spectral characteristics of 3-min oscillations above the umbra and 
tendency of their spatial distribution we will start from the generally accepted idea 
that these oscillations are caused by magnetoacoustic gravity (MAG) waves, that is the 
acoust ic-like slow mode propagation along magnetic field lines. It has been long known 



(since 



Lamb 



1908! ) that due to the presence of gravitational field compressible acoustic 
waves can propagate upwards through the temperature-minimum only at frequencies above 
the cut-off frequency: 

where g is gravitational acceleration and 

Cs = VlRT/fi (3) 
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is the speed of sound, T is the temperature, 7 is the ratio of specific heats, fi is the mean 
molecular weight, and R is the gas constant. 



Further, 



Bel fc LeroyI (11977 ) showed that for MAG waves in the solar atmosphere 



the cut-off frequency will be aflFected by the magnetic field when it is not parallel to the 
gravitational field. They predicted that in regions of low-/3 plasma (/3 << 1), which is the 
case above the sunspot umbra, the effective gravity on a particular magnetic field line would 
be decreased by the cosine of the inclination angle 9 of that field line with respect to the 
direction of gravity (the solar normal), that is 

g = QqX cos9, (4) 

where go = 274 m s~^. Therefore, cut-off frequency 

fc oc cose/Vf. (5) 

Based on the above theoretical predictions one may expect that the change of the observed 
cut-off frequency across the umbra may result either from variation of temperature in the 
temperature-minimum layer, T^^ri, becoming smaller in the center of the umbra or the 
magnetic field inclination angle 9 becoming large enough near the umbra/penumbra border. 
The upper limit of fc is determined by T^^^ on the path of the MAG wave propagating 
along the vertical field line. 



According to Figure [5] the value of cut-off fre quency in the center o: 



the umbra, in 



exactly vertical field, is fc = 5.7 mHz. If, following iFleck fc Schmitd ( ll99ll ). we take /x = 1.3 
and 7 = 5/3 for an ideal neutral gas, using Eq. [2] we obtain T^in ^ 3800 K. The uncertainty 
of the this estimate is mainly determined by the uncertainty of the mean molecular weight 
jn: change in /x by 0.1 gives the change in temperature 300 K. Moreover, the mean molecular 
weight in turn depends on temperature. 

Near the umbra's boundary observed /c = 5 mHz. Using proportionality [5] under the 
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assumption of T^in = const gives the inclination angle 6 ^ 30° at the peripheral parts 
of the umbra. This value is consistent with the results of previous study of the radial 
variation of field inclinatio n on the level of the temperature minimum (e.g. Fig. 3 in 



Mcintosh fc Jefferies 



20061). 



On the other hand, it is reasonable to assume that the minimum temperature grows 
from the middle to the peripheral parts of the umbra. Under the assumption of vertical 
field through all the umbra the observed value fc = 5 mHz Rives Tm.j^, ^ 4900 K. This 
value is too large even for the quiet-Sun atmosphere ( Avrett fc LoeseL.2008i ). Therefore, we 
think that the temperature increase has a smaller effect on the wave frequency distribution 
across the umbra than the magnetic field inclination, though both effects modify the cut-off 
frequency in a similar way. 

The gradual broadening of the area occupied by the oscillations with the height of 
the solar atmosphere agrees well with the above interpretation. Roughly speaking, those 
field lines inclined 30° to the solar normal must form a wider circle in the transition region 
than at the temperature-minimum level. An estimate using the simplifying assumption 
that the magnetic field remains linear between the two formation heights (there is no field 
curvature), gives 1-2^' of the circle expansion. This estimate well agrees with the observed 
value Ad shown in Section I 



requency modificatio n in t 



ne 



Mcintosh fc JefferiesI ( 120061 ). 



Observational demonstration of the acoustic cut-off 
penumbra area using TRACE observations presented by 
Results obtained in this work can be considered as a demonstration of the same effect in 
the umbral area. Nevertheless, magnetic field extrapolation would be useful to confirm it. 

The other possible explanation of ring-like frequency-sp atial structure can be done 



in th e frame of the chromospheric acoustic resonator model (I Zhugzhda 



20081 : 



Botha et al. 



20111). Indeed, the increase in the angle of the guiding magnetic field from vertical, the size 
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of the resonator grows and hence the resonant frequency decreases. 

5. Conclusion 

Three-min oscillations above a sunspot have been studied both at short (UV/EUV) 
and long (microwave) wavelengths. Intensity variations with / = 5 — 9 mHz are found 
in all analyzed UV/EUV bandpasses as well as in microwaves (1.76 cm) with the similar 
spectral characteristics. The largest modulation depth (20%) is found in the transition 
region (304 A channel) and it decreases down to the photosphere and up to the corona. 
Time delays between oscillations at different wavelength above the umbra's center increase 
with the expected formation height of the corresponding emission (except for 211 A). This 
confirms that observed intensity variations of both UV/EUV and microwave emission result 
from waves propagating upwards. Measured delays are used to estimate average velocities 
of wave propogation at different layer of solar atmosphere and the altitude of microwave 
emission generation. 

High spatial, temporal and spectral resolution of SDO/AIA allowed us to trace the 
variations of the cut-off frequency and spectrum of oscillations across the umbra. We found 
that in the range 5-9 mHz high-frequency oscillations are more pronounced closer to the 
umbra's center, while the low-frequencies concentrate to the peripheral parts. Accordingly, 
the cut-off frequency decreases from 5.7 mHz to 5.0 mHz. We interpreted this discovery 
as a manifestation of variation of the magnetic field inclination across the umbra at the 
level of temperature-minimum. This interpretation has been used for diagnostics of sunspot 
atmosphere on this level. 

The similar spectral characteristics of 3-min oscillations were found through all the 
temperature bandpasses. This suggests that the strong magnetic field of a sunspot works as 
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a waveguide for the acoustic waves propagating from the photosphere level and eventually 
reaching 1 MK corona. 
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Fig. 1. — (a) SDO/HMI magnetogram of the sunspot region in NOAA AR 11131 and images 
of sunspot taken (b) by NoRH and (c-h) by SDO/AIA on 2010 December, 8 at 04:00 UT. 
Orange ellipse at top-right of plot (b) indicates the half-maximum NoRH beam size. White 
cross and asterisk on plot (c) indicate two pixels inside the umbra for later reference. 
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Fig. 2. — (a) The time profiles of intensity above the center of the umbra from six AIA 
bandpasses and Tsmax at 1.76 cm of 22 min duration. Emission modulation depth at (b) 
1700 A (solid), 1600 A (dash), 304 A (dash-dot-dot), and 1.76 cm (dash-dot) and (c) at 
171 A (thick), 193 A (dots), and 211 A (thin). For the color version of the figure, different 
wavelength are shown by colors, corresponding to the key to the right. 
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Fig. 3. — (a) Cross correlation coefficients of long time series of pulsations above the umbra's 
center of each analyzed wavelength with respect to 304 A as a function of the lag. The same 
colors and line types as in fig. [2] are used to indicate different wavelength, (b) Time lags 
obtained from cross correlation for each wavelength with the corresponding uncertainties 
shown by error bars. 
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Fig. 4. — Emission modulation depth of signals at 304 A (red) and 193 A (black) from the 
center of the umbra filtered in narrow (1 mHz) frequency bands with the central frequencies 
at 5 mHz, 6 mHz, 7 mHz and 8 mHz (± 0.5 mHz). Signals at two frequencies are overplotted 
with a time shift 6t to match each other. 
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Fig. 5. — Fourier power spectra of the intensity variation above the center of the umbra for 
all analyzed wavelengths. 
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Fig. 6. — Fourier power spectra of the intensity variation near to the umbra/penumbra 
border from four AIA bandpasses. 



-29- 




Fig. 7. — Spatial distribution of normalized Fourier power in four frequency bands with the 
central frequencies at 5 mHz, 5.5 mHz, 6 mHz, and 7 mHz (±0.2 mHz in each band). Power 
grows with brightness. The umbra-penumbra boundary is shown by contour line, white on 
panels (a-b) and black on panels (c-f). 



